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Abstract 

We propose some useful estimates for the pointwise error estimates 
of the streamhne diffusion finite element method (SDFEM) on Shishkin 
meshes, when SDFEM is applied for problems of characteristic layers. 

1 Problem 

We consider the singularly perturbed boundary value problem 

,^ -eAu + bux + cu = f inf7 = (0,l)^, 

^ ' u = on dn 

where 6, c > are constants and 6 > /? on f2 with a positive constant /3. It is 
assumed that / is sufficiently smooth. Here < e ^ 1 is a small perturbation 
parameter whose presence gives rise to an exponential layer of width 0(e) near 
the outflow boundary at x = 1 and to two characteristic (or parabolic) layers of 
width 0{^/e) near the characteristic boundaries at y = and y = 1. 

2 The SDFEM on Shishkin meshes 

In this Section we describe our mesh, our finite element method and the as- 
sumptions of our analysis. 

2.1 The regularity result 

As mentioned before the solution u of (jl.ip possesses an exponential layer at 
X — 1 and two characteristic layers at y — and y = 1. For our later analysis we 
shall suppose that u can be split into a regular solution component and various 
layer parts: 
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Assumption 2.1. The solution u of can be decomposed as 
(2.1a) u = S + E1+E2 + E12. 

wliere for all a; = (x, j/) G O and for < i + j < 3, the regular part satisfies 



(2.1b) 



< c, 



while for the layer terms and < i + j < 3, the following bounds hold true: 



(2.1c) 



■ix,y) 



(2.1d) 

and 

(2.1e) 



d'+^E2 



dx^dyi 



ix,y) 



For constant coefBcients Kellogg and Stynes [H [S] give sufficient compatibil- 
ity conditions on the data that ensure the existence of (|2.1ap ~ (|2.1ep . 



2.2 Shishkin meshes 

When discretizing (jl.ip , we use a piecewise uniform mesh — a so-called Shishkin 
mesh — with iV mesh intervals in both x— and y— direction which condenses in 
the layer regions. For this purpose we define the two mesh transition parameters 

Ax := min I -, p— In j- and Aj, := min | -, /3\/£ In 

In this paper, we define p — 2.5. 

Assumption 2.2. We assume in our analysis that e < . Furthermore we 
assume that A^^ = p£l3~^ InN and Aj^ = p^/e\nN as otherwise N^^ is exponen- 
tially small compared with e. 

The domain 51 is dissected into four(six) parts as 51 = fi^ U 51i U U 51i2, 
where 



[0,1- A,] X [Aj,,l- A, 



fly 



fli [1 - Xx, 1] X [Aj,, 1 - Aj^] , fli 



[0,1 -A,] X ([0,Aj,]U [1- Ay,l]), 
= [1-A,,1] X ([0,A,]U[1-Ay,l]). 



Remark 2.1. The mesh transition parameters have been chosen such that 
the boundary layer function E which can be any of Ei , E2 and £'12 satisfies 
\E\ < CN-P on VI,. 
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We introduce the set of mesh points {{xi,yj) £ fl : i, j — 0, ■ ■ ■ , N} defined 

by 

2i(l-A^)/7V, for i = 0, ••• , iV/2, 

1 - 2{N - i)X^/N, for i = N/2 + 1, ■■■ , N 

and 

r ^j\/N, for J = 0, • • • , N/3, 

Vj^l {3j/N^l)^3{2j^N)\y/N, for j = 7V/3 + l, ••• , 2iV/3, 

[ 1 - 3(7V - j)Xy/N, for j = 2N/3 + 1,---,N. 

By drawing Unes through these mesh points parallel to the x-axis and y-axis the 
domain fl is partitioned into rectangles. This triangulation is denoted by fl^ . 
If D is a mesh subdomain of il, we write for the triangulation of D. The 
mesh sizes h^^r = Xi — Xi-i and hy -^ = Vj — Vj-i satisfy 



^^^^IV/^' iovi^l,--- ,N/2, 
hx := for i = N/2 + !,■■■ ,N, 



N/2' 

- 2X, 
N/3 



1 — 2A 

Hy~ for j = 7V/3+l, ••• , 2iV/3, 



hy '■— ,^^„ , otherwise 



iV/3' 



and 



N-^ < Hx,Hy < 3N-\ 
CieN-^ \uN <K< C2eN~^ \nN, 
Ciy/eN-^ \nN < hy < C2\/eN-^ In TV. 

The above properties are essential when inverse inequalities are applied in our 
later analysis. 

For the mesh elements we shall use some notations: Tij = [xi-i, Xi] x [j/j-i, yj\ 
for a specific element, t for a generic mesh rectangle (see Fig[T]) and 



{x^^i+x.i)/2, yr ^ ivj-i + yj)/2 if t = [x,-i,Xi] x [yj^i,yj] 



2.3 The streamline diffusion finite element method 

The weak formulation of the problem is: Find u £ -ffg (fi) such that 

(2.2) e(Vu, Vf) + {bux + cu,v) = (/,«), Vw £ H^{^). 



Note that the variational formulation (j2.2p has a unique solution by means of 
the Lax-Milgram Lemma. 
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l3 



l1 



Fig.l: Geometry of the element r 

On the above Shishkin mesh we define a finite element space 

:= {v^ e CiCl) : u^lao = and v^\r is bilinear, Vr G Vl^}. 

Then we can state the standard Galerkin discretisation of (|2.2p is: Find U € 
such that 

(2.3) e(VC/, Vv^) + {bU^ + cU, v^) = (/, w^), Vw^ e F^. 

The SDFEM adds weighted residuals to the standard Galerkin finite element 
method : Find U eV^ such that 

(2.4) B{U,v'') = {f^v'^+Sbv^), V^^ e F^, 
where 

(2.5) B(f/, w^) £(V?7, Vw^) + (fef/^ + cU, v^) + {bU^ + cU, Sbv^). 

The term {—eAU,6bv^) is neglected in our case. S — 6{x) is a user-chosen 
parameter (see [Hl^). In this paper, we set 

S{x) <^ 

L 0, otherwise 

where C* is chosen so that < S{x) < l/cfoi x & flgU fl2 (see [3 §111 3.2.1]). 
Finally, we define a special energy norm associated with _B(-, •): 

(2.6) \\\U\\\^ {{e + b^S)U,,Ucc) + e{Uy,Uy) + {cU,U). 

For any subdomain _D of fl, let Bd{-,-), {■,-)d and ||| • |||£) mean that the 
integrations in (|2.5p and (|2.6I) are restricted to D. Wc denote by || • ||u the 
norm in L^{D), i.e., 

||u||2^ = {v,v)d for aU w e L^{D). 
If D ^ n then we drop from the notation. 
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3 Interpolation error estimates 

We start our analysis by quoting some previous results. In the following analysis, 
we shall frequently use the bilinear interpolation of a given function g. 

Lemma 3.1. Let t G Vl^ and p G [l,oo]. Assume that g lies in C^(t). Then 

IKff - g^)x\\LP(T) < C {hl.^^\\gxxx\\LP(T) + hx,Thy^r\\gxxy\\LP{T) + hl ,r\\9xyy\\LP{T)) 

+ Chx^T\\9xx\\LP{T)j 

IKff - g^)y\\LP{T) < C {hl^^\\gxxy\\LP{r) + hx^rhy^rhxyyW Lp{t) + ^l^rhyyvW Lp{t)) 
+ Chy^r\\gyy\\LP{T)- 

Proof. See T, Theorem 4]. □ 

Lemma 3.2. Assume that u satisfies Assumption 2.1. Then on our Shishkin 
mesh 



\u-U^\\L'^(r) < 



CiV"^ In^ otherwise. 



Proof. See % Theorem 3]. □ 

Lemma 3.3. For any function g G C^{t) and any w G , we have the 
identities 

(3.1a) j {g - g')xWa;dxdy J g^yyj^{y) - ^(y - yr)wxy^ dxdy, 



(3.1b) 
and 
(3.1c) 
where 



^ {g- g')yWydxdy = J gxxyFr{x) (^Wy - ^(x - Xr)wxy^ dxdy 



J (9~ g')xwdxdy = J R{g, 'w)dxdy + (/ ^ / ) 5^^"^^^ 



Fr{x)^\{(x~Xrf~hlJA) and My) ^\{{y - Vrf ~ hlJA) 

and 

1 Kt 
R{g,w) ^-Fr{x){x - Xr)gxxxWx ^gxxxW + Jr{y)gxyy 



2 2 

w- {x- Xr)wx - - yT)wy + 3(2; ~ ^T){y - yr)Wxy 



Proof. See [BJ [7] or [H the Appendix] for details. □ 
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Lemma 3.4. Let U be the solution of (|2.4I) on our Shishkin mesh. Then 

\\\u^ -U\\\< C{N-^ N + e^/^N-^ In^/^ 
Proof. See 2, Theorem 5]. □ 

Throughout the remaining analysis we shall make frequent use of the follow- 
ing inverse estimates. Let x be a polynomial on the mesh rectangle r. Then 



(3.2) 
(3.3) 
(3.4) 



\\Xx\\lv{t) < Ch^^^\\x\\LP{T), WXyhpir) < C/l^^J. || x|| LP(r) , 



\xixi,y)\dy < Ch^]^\\x\\LHr,,), 



Vj-i 



\\x\\LHr}<CiK^rhy,ry^''-'/P\\x\\LPir) for P,q e[l,Oo] 



Lemma 3.5. Let Assumption 2.1 hold true. Then there exists a constant C 
such that the following interpolation error estimates hold true 



(3.5a) 
(3.5b) 
(3.5c) 
(3.5d) 
(3.5e) 



{E,-Ei)y\\L^^n,)<CN-'', 
Viu-u')\\L^in.)<CN-\ 



Ei-Ei\\L^(^n,)<CN-P, 

{Ei2 - Ei^UlL^^n,) < Ce-'N-P, 

iE,2~Ei2)yh^^n.)<Ce-'/'N-P, 

{E,2-Ei2U\LHn,)<Ce'/^N-P, 

Viu-u')\\Liin,)<CN-HnN 



where the function E can be any one of Ei , E2 or £^12 . 

Proof. At first, we will prove the second inequality of (|3.5dp . The proof of 
p.5ap - (|3.5c|) and the first inequality of p.5dp is similar. 

Each bilinear basis function (f)^'^[x,y) satisfies <jf'^{x^y) — (t)i{x)(j>' (y) where 
(jji and (f)^ are piecewise linear basis functions. If a; e [xi^i, Xi], 

I / \ Xi X 1 / \ X Xi — i 
(pt^l[x) = — , 01 (X) = . 

The functions (j)^~^{y) and (/>■' (y) are defined similarly in ?/-,■]. We define 

z{x, y) := Ei2{x, y) and w{x, y) := e-''(i^^)/^(e^^^/^ + e-^^-y^^). 
A direct calculation and (2.1e) give 

(3.6) ||(-Bi2).||lmo.) < Ce-^\\w{x,y)\\L^(n,) < Ce^'^N^P. 
For {El2)x^ we have 

{Ei2)x\r., {z%\,., = £.^.i_£^M^^.-i(2;) + £M_£^^.(y) 

(3.7) = h-]^ i^^^iy) z^ix, yj-i)dx + LfP (y) z^{x, yj)dx^ 
= K}t [Ixi, ^^(^' yj-i)<^x + if^ (y) J^^_^ z^y{x, t)dxdtj . 
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Similarly, we have 

(3.8) (zO.lr,, = K}, (j^l^ z,{x, y,)dx - v^-\y) JJ^^^ 

By direct calculations, we can obtain 



Zxydxdt 



(3.9) 



w{x,y) > w{x,yj) > if [yj-i,yj] C [0,Xy], 
w{x,y) > w{x,yj-i) > if [yj^i,yj\ c [1 - Ay, 1] 



for y G [yj_i,yj\. Then, for any Tij G Vty and (x, y) G r^, from (2.1e) we have 

(3.10) \z^y{x,y)\<Ce-'"^w{x,y) 
and 

(3.11) \zx{x,yj)\ < Ce^'^w{x,y) if C [0,Xy] 

(3.12) |z^(a:,yj_i)| < Ce"iw(x,2/) if [y^.i, y^] C [1 - Aj^, 1] 

where we have used ^^9^. Combining ({371) . ([3J0)) . ((3l^ or jSH), ([XTO)) . (jXTTI) 
and considering < (j>'~^{'y), 4^ {y) < 1, we obtain 



\\{Ei2)A\L^i.n,) = |(^').|da:dy 



< C 



e "'^w(x, y)da; + 



z ^/'^w{x,t)Ax(lt I dxdy 



Then, we have 

< C{e-' + e'/'N-' \nN ■ e-'/')\\w{x, y)\\mn,) 

From p. lap , we have 

u-u' = {S^S') + {El - Ei) + {E2 - E^) + {Ei2 - Ei^). 

For the proof of the first inequality of (j3.5ep . from the standard interpolation 
theory and (|2.1bp . we have 



(3.13) 



l|V(5 - 5^)|Ui(o.) <CiV-^ E 

i+j=2 



dx'^dyi 



< CN-\ 



(|2.1cp and the inverse estimates p.3p give 

(3.14) ||V(i?i - i?[)|lLi(J7.) < ||V£;i|Ui(o.) + WVEiWmn.) 

< l|V£;i||Li(n.) + CN\\EiUiin^^ < CN'-p 
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where we have used < CN p for {x,y) € Vts- SnTiilarly, we have 

(3.15) \\V{E2 - Ei)\\Li^n^) + \\V{Ei2 ~ S(2)IUMf^.) < CN^-p. 

Combinmg (jXTl) . ([5J5)) . we are done. 

For the estimate of the second inequahty of (I3.5ep . on the one hand, we 
apply Lemma ED to V{S - S^), V{Ei - E{), V{E2 - E^) and {E12 - £^(2)^ 
and on the other hand, apply the similar analytic techniques as for the second 
inequality of (j3.5dp to {E12 - E{2)x- □ 

4 The discrete Green's function 

In this section, we will introduce the discrete Green's function and derive some 
estimates of it . 

Let X* = {x*,y*) be a mesh node in The discrete Green's function 
G G V'^ associated with x* is defined by 

(4.1) B{v^,G)^v^{x*),yv^ eV^. 
We introduce 

(4.2) = kN-^ In TV, ay = kN-^''^. 

The constants fc > 0, sufficiently large and independent of N and e, are chosen 
according to the derivation of Theorem 14.11 

Theorem 4.1. For a;* e fi^ U ili we have 

\\\G\\? < 8\\\G\\\l < CNlnN. 

Proof. See [HI Theorem 4.1] □ 

For pointwise bounds on G and its first-order derivatives, we define a sub- 
domain of fl as 

flo ■■= {x e fl : X ~ X* < K<Tx\nN and \y~y* \ < KaylnN}. 

The constant K > will be chosen later. 

We extend ilo to the smallest mesh domain H'q — ilQ{x*), i.e., 

% ^ U{t e : meas(r2o n r) ^ 0}. 

Note that meas(r2Q) < Gay \nN. 

Theorem 4.2. Assume that a^ — kN~^ \nN and ay — kN^^^^, where k > is 
sufficiently large and independent of e and N. Let a;* e fig U ili, then for each 
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nonnegative integer v, there exists a positive constant C = C{v) and K = K{v) 
such that 

\\G\\w^-°°(n,\n'g) < CN'"" , 
e|G|Tyi.~(ni\n^) + ||G||L~(ai\f2^) < CN^"", 
e|G|vKi.~(ni2\n;,) + l|G||L~(ai2\f2;,) < Ce^^'^N'^ 

and 

+ e^/%GyU^(n.\n'„) + ei/'l|G|U^(a,\n;,) < CN-\ 
Proof. The proof is similar to the one of iSj Theorem 4.2]. □ 

5 Maximum-Norm error estimates 

In this section we shah derive bounds on \ {u — U){x)\ for x lying in the various 
subregions of 51. 

Taking = U - \n (|lT|) yields 

([/ - u){x*) = (f/ - u'){x*) = B{U - u^, G) 

where x* is a mesh node. Let e denote the interpolation error, i.e., 

e{x) := (u — u^){x). 

Then from ([221) , (1231) and we have 

(5.1) {U - u){x*) = -£(Au, SbGo,) + S(e, G). 

The various terms on the right-hand side are bounded separately. 
For the following analysis, we derive some useful local estimates. 

Lemma 5.1. If Assumption 2.1 hold true, then there exists a constant G such 
that 

(5.2a) WiE ~ E'U^,^^^^^,^ < CN-Pay\nN, 

(5.2b) \\{{E^ + E,2) - {Ei + Ei,))y\\^,^,,^^,,,^ < CN-Pay\nN, 

(5.2c) \\A{Ei+Ei2)\\LHn,nn'„) < Ce-^N-PaylnN 

where ay as in (|4.2p and the function E can be any one of Ei, E2 or Ei2. 

Proof. The proof of (|5.2ap and (j5.2bl) is similar to [101 Lemma 4.1]. Inequality 
(j5.2cp can be deduced directly by Assumption 2.1 and the definition of JIq. □ 

Theorem 5.1. Assume that ax = kN^^hiN, ay = kN^'^/^ and e < N^'^ . 
Then, for x* Cz fls U ili, we have 

|S(e,G)l < C{N-^/^ + e^^^N-^){lii'^ N) ■ \\\G\\\. 
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Proof. For the following analysis, we define 

Ci {{n, u fii) n n'o) u {^2 u ^12) 

and modify the bilinear form by mean of integration by parts and the decom- 
position (|2.1al) as follows: 

B{e, G) =((£ + b''S)e^, G^) + e{ey,Gy) + {b{S ~ S')^,G) 

+ {b{E2 ~ eI)x, G)n2uai2 + {HE2 - E^), Gux) g^n.uni) 

— {b{E - E^),Gx)n^uni - {b{Ei — E{), Gx)n2UOi2 

- {b{E^2 - Ei^), G^)n,^jn,, + (ce, 5bG^) + (ce, G) 

where E = Ei + E2 + -E12, is the x-axis coordinate of the outward normal 
vector of 9(f2s U ili). Note that (fe(i?2 — E^), Gnx)d(n^uni) — because G — 
on d^l. 

The discussion of i3(e, G) will be separated into three parts. In (a), we will 
analyze ((e + b'^S)ea;,Gx) and e{ey,Gy). In (b), {b{S - S^)^,G) and {b{E2 ~ 
-^2)2:1 G')o2usii2 ■'^ill be discussed. In (c), we will analyze the residual terms of 
B{e,G). 

(a) In this part, based on the boundary layer behavior of u, we discuss u — 
by means of the following decomposition 

u-u^ ^{S-S^) + [El - Ei) + {E2 - E^^) + {Ei2 - Ei^). 

According to Lemma |3.3[ we have 

{{S - S^)x, Gx)^ = j SxyyJriy) {gx - ^(y - yT)GxT^ dxAy. 

Then 

I ((S* — S Gx) J\ < C/l^^^llSajj^j^llT-IIGajllT- 

where we have used the inverse inequalities (13. 2p . Thus, 

(5.3) \{{e + bH){S-S\,Gx)^\ 

< CN-me + bH)^/^Sxyy\\^,\\{e + b^S)^/^G,\\n 

< CN-^/^{aylnNy/^\\\G\\\ < CN-^^^^{ln^/^ N)\\\G\\\ 

where we have used meas(f2) < Cay In + Ce^/^ In < Cay In A^. 
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Similarly, we have 



{{e + b^S){E,-ElUG,)^^^^,\<CN''/\\n'/'N)\\\G\\\, 

{is + bH)iE2-EiUG,)^^^^^^\<Ce'/^N-'^/\\n'N)\\\G\\\, 
{ie + bH)iE,,-Ei,),,G,)^^^\<Cs'/^N~'{\n'N)\\\G\\l 
\{{S~S%,Gy)^\<Ce'/'N-'/\\n'^'N)\\\G\\\, 



(5.4) 
(5.5) 
(5.6) 
(5.7) 
(5.8) 
(5.9) 
(5.10) 
(5.11) 
(5.12) 

Furthermore Lemma 15.11 and the inverse inequality p.4p imply 

(5.13) |((e + b^5){E - E').,,G.,)n^nni, \ 

< CN ^\\{E — E^)x\\L^{n^r\%) ■ \\Gx\\L'=<'{n^nn'g) 

< CN-'N-Pay{\nN) ■ ^||G,||o.noi 
<C7V-''(ln7V)|||G|||. 

Similar argument shows 

(5.14) e\{{E,-Ei)y,Gy)^^^^,\<Ce'/'N'/'-'^i\nN)\\\G\\l 

(5.15) £|((i?i2-£;(2)„G,),,^^^,| <Cei/2iVi/2-P(iniV)|||G|||, 

(5.16) {{e + b'6)iE^2~Ei^)x,Gx)^^ < Ce'/^N'^'-P{ln-'/' N)\\\G\\\ 



{{E2 - E2)y,Gy)^^^^^^,^ 

{(.E~E%,Gy)^^^^^^,\ < Ce'/'^N-^W N)\\\G\\\, 
{{E-E')y,Gy)^^^\<Ce'/^N-\\n'N)\\\G\\\, 
{{E2-Ei)y,Gy)^. <Ce'/^N 



where we have used p.5dp in the last inequality. 
Lemma [3.51 and Holder inequalities give 



(5.17) 



{ie + b'S)iE2-Ei)x,Gx)^^^^, 

<£\\{E2 - E2)x\\L°°{ninn'„)\\Gx\\Li{ninn'„) 
<CeN-f{e<jy\n^ Ny^^\\Gx\\n,nn', 
<GeiV-(i/4+'')(ln^)|||G|||. 
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Similarly, we have 

(5.18) |((£ + fe^5)(£;i2-i?ig.,G.)^^^^,| <C7iV-(i/^+'')(ln7V)|||G|||, 

(5.19) e\{iE,-Ei)y,Gy)^^\<Cs'/^N-Pi\n'/'N)\\\G\\\, 

(5.20) e\{iE,,-Ei,)y,Gy)^^\<Ce'/^N-P{\n'/'N)\\\G\\\. 
In view of Theorem 14.21 with v — 2 and Lemma 13.51 we see 



(5.21) I {{e + b^5){u - ■u^)^;, Gx)(n,uni)\n^ + £((« - u^)y, Gy)(^n,uni)\n'g \ 

<CN-^. 
(b) We see from Lemma [3731 that 

((5 - S'h,G) f ^(^' + Y,^{I - f) S^^Gdy 

ren-'^ ren ^-^'^ -'(4/ 



where i?(-, •) as in Lemma 13.31 Based on our Shishkin mesh and the properties 
of the discrete Green function G, we decompose the first term as follows: 

J2 j R{S, G)dxdy = ( ^ + ) j Ri{S,G)+R2 {S, G)dxdy 

where 



1 

R,{S,G) = -FAx){x - ..)5_G. - ^^_G 



and 



2 2 

i?2('5', G) = JT{y)Sxyy{G~ {x - Xr)Gx ~ gCj^ ~ yr)Gy + 3 (^^ " a:^r)(y " yT)Gxy). 



Firstly, from Assumption 12.11 and the definition of Q, we have 
(5.22) ^^l|F,(x)(a;-Xr)5 

XXX 

Gx\dxdy 

Ten 



Ten 

< CH^\\Gx\\L^((n^n%)un2) + Ghl\\Gx\\Lmninni-,)uni2) 

< GN-^^/\ln^/^ N)\\\G\\\ 



and 



(5-23) E / ^\Sxx.G\dxdy < CN-^/^{ln'^' N)\\\G\\\. 



Ten 
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Applying inverse inequalities p.2p to the last part of i?2(5', G), we obtain 

(5-24) Y.Jj^2{S,G)\dxdy <Cj2hlr\\S.yy\\r\\G\\r 

< GN-"^ ■ {uylnNY'^ 

< GN-^/\\n^/^ N)\\\G\\\ 

From Theorem 14.21 with v — 1, we have 



(5.25) J2 I \Ri{S,G)\Ax<ly 

< C ^ [hl^^\\Sxxx\\L^T)\\Gx\\L-°(T)+hl.^^\\Sxxx\\L^{T)\\G\\L-°{T)) 

rG(i2sUOi)\n^ 

< CHl\\Gx\\L^(n,\n'a) + Ghl.\\Gx\\L^(ni\Q.'^) + C'^~^llG|lL=-((a,ur!i)\o^) 

< CN-^ 

and 

(5.26) [\R2iS,G)\<G J2 H^\\Sxyy\\LHrm\L^ir) 

Te(o.,uf2i)\f2^ Te(i2aUOi)\o^ 

< GN-^ 



where we have used inverse inequalities p.2p . 

Secondly, we set L := {{1 — , y) : < y < 1}. Then we have 

T.^(f - f) ^-Gdy = 1 ~ hi) f SxxGdy 

<Ci/J f l^-^l < CH^ I ^ / dy+ Y f ] 

/6-L ' \lGLn% ' lGL\n'o^ J 

^GH^Jl + ll). 
The estimate of I is straightforward: 

1= E / \s..G\(i-xx,y)dy< r r 

< \\Sxxx\\l°°{Dl)\\G\\l^Dl) + \\Sxx\\l°°{Di^)\\Gx\\l^Di,) 

<G{eay In' N)'/^{\\G\\n. + \\Gx\\D.) 

< GN'^/HnN- {e'/^G\\D,+e'^'\\Gx\\Dj 

< GN-^/\lnN)\\\G\\\ 



{SxxxG + SxxGx)dx 



dy 
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where {1 - A^} x [y[,y'2] = Lnfi'o and Dl := [1 - A^, 1] x [y[,y'2]. 
From Theorem 14.21 with u = 1, we have 



\SxxG\ dy 



n= \S..G\dy+ 

< C||G|U^(oAn;,) + Cei/2(in7V)||G|Uoo(o,\o;,) 

< CN-\ 

Considering the estimates for I and II, we obtain 



(5.27) 
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Ten 



1-2 <J I4 



SxxGdy 



< CN-^^\lnN)\\\G\\\. 



The estimates for {{E2 — E{)x,G)n2uni2 are the same as {{S — S^)x,G). 
Thus we have 

(5.28) \b{{E2-E^)x,Gh2un,2\<Ce'/^N-^{\n^N)\\\G\\\. 
(c) From Lemma 13.51 we have 

(5.29) \b{E~ E', Gx)^^^^,\ <C\\E- i?'||L-a2.m2;,)l|G.|Ui(i2.nnf,) 

< CN-P ■ {ay\nNY/^Gx\\n^nn'„ 

< CN^^^-P{ln^/^ N)\\\G\\\. 

Similarly, we have 

(5.30) \b{E,-Ei,Gx),,^\<Gs'/^N'/'-P(\n'/'N)\\\G\\l 

(5.31) \b{E,2'Ei2,Gx)^^\<Ce'/'N'^^-f'{\n'/^N)\\\G\\\. 
In view of Lemma 13. 2[ we obtain 



(5.32) ^{^ ^ ^Gx)^^^^,^ < G\\E ~ E'\\L^(^Q^r]n'g)\\Gx\\L^ninn'„) 

< GN-^{\n^N) ■ {eayln^ Ny/^\\Gx\\n,nn', 

< GN-^/\\n^ N)\\\G\\\. 

Similar argument shows 

(5.33) b{Ei- Ei,Gx)^^^^ < Ge'^/'^N'^^{\n^ N) 

(5.34) \b {E,2 - Ei2,Gx),,^^ \ < Ce^/^N-^ln^ N) ■ 

(5.35) \c{u~-u\5Gx)\ < C7V-^/^(ln^7V)|||G|||, 

(5.36) \c{u — u , 
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Theorem 14.21 with v ^ 1 and Lemma \S72\ yield 

(5.37) \b{E-E',G,)^^^^, 

< C\\E - £;^||i^(aj)meas(f^i \ ^'o)\\G,:\\L'^(^m\%) 

< CN-^ In^ N{e\nN)\\G4L^^n,\n',) < CN~^ 

and 
(5.38) 
(5.39) 

CoUecting (j5.3p - (|5.39l) . we are done. 



{u~u',G)^n^un,)\ni, \ <GN-\ 



□ 



Theorem 5.2. Assume that u satisfies Assumvtion \2.1\ and e < N ^ . Let 

= kN-^^ In TV and ay = kN~^/^. Then for x* efl^U fli 

1. if ^I'o C U ill, we have 

\ieAu,SbG^)\ < G{N-P + eN-^/^)i\nN)\\\G\\\. 

2. if ftgU fti, we have 

I (sAu, 6bG,)\ < C{N-P + £1/4(5,) (In iV) 1 1 |G| 1 1 . 

Proof. We set E = Ei + E2 + £'12 and define Tg^x '■— H fii and Ty^^y '■= 
U,2 n r2i2. At the beginning, integration by parts and the definition of 5 yield 

{AS,5bGx) =iAS,SbGx)n^un, 

^{AS,S,bG)r^,^ + {AS,SybG)r^^^^ ~ ((A5),,56G) 



and 



iAE2,SbG,) = {AE2,SsbG)r^,_^ + iAE2,SybG)r - {{AE2)^,5bG). 



Thus, 



{eAu, SbGx) =eSsiA{S + E2), 6G)r,,, + sSyiA{S + £;2), &G)r„,,, 
- £((A(5 + E2)):c,SbG) + e {A{Ex + Ei2),SbGx) 



The terms on the right-hand side are analyzed separately. 
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From Assumption 12. 1[ we have 
(5.40) eS, \{AS,bG)r^,J < CeN-^\\AS\\L^^r^.J\G\\mr^^^) 



<CeN-'l / |G(l-A.,y)|dy+||G|Ui(r_\o^) 
<CeN-^lf [ \G,{x,y)\dxdy + \\G\\L^^a^\n'j] 

< GeN-^/^{liiN)\\\G\\\ + GeN-^ 
where we have used Theorem 14 . 2 1 with v = 1. Similarly, we have 

(5.41) ed,\{AE2,bG)r^J < GN-^'+p\\n'/^ N)\\\G\\\, 

(5.42) sSy\iAiS + E2),bG)r,,J < Ge'/%i\nN)\\\G\\\. 

Considering (j2.1bp and meas(r2 n fig) < Gay IniV, we obtain 

(5.43) \s{{AS),,SbG)^^,,,\ < GeN-'\\iASUL^^nnn'J\G\\mnnK) 

<GeN-\ay \nN)^/^G\\nnn', 
< GeN-^^\lii^/^ N)\\\G\\\. 

Assumption 12.11 and the inverse inequality p.4p yield 



(5.44) \eiiAE2).,6G)n,nn'„\ < Ge6s\\iAE2)o.\\ )\\G\\l 

< GeSs ■ e-^I^N-" ■ {H^Hy)-^l^\\G\\a^r^n', 
<e^l^N'P ■\\\G\\\. 

In view of (j2.1dp . we get 

(5.45) \e8y[[AE2UG)^,r,n',\<Ce8y\\[AE2\A\i-i.\\G\W 

< Ge^l^K 



Lemma 15.11 Assumption 12.11 and the inverse inequality p.4p yield 

(5.46) |£(A(i?i+i?i2),<56G,)^,,nOi 
<Ge5,\\A{Ei+Ei2)\\ 

<Ge5,-e'^N-fay\^N ■{H,Ily)-^l''\\G^\\a^r^a', 
<GN'P{\^N)\\\G\\\ 
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and 

(5.47) 



n2nsif, 



<CeSy\\A{Ei + Ei2)\\Li{n2n%)\\Gx\\L'^{n2n%) 
<Ce6y ■ e-'/^N-P{\nN) ■ {H^hy)-'^^\\G,\\n2n% 
<Ce^/'^N^/^-P{ln^/^ N)\\\G\\\. 

Considering (|2.1b[) — (|2.1cp and Theorem 14 . 2 1 with w = 1, we obtain 



(5.48) 



■((A5),,<5feG)^^^, < CeS\\{ASUL^(n\n'JG\\L^^n\n'j 



Similar argument shows 
(5.49) 



\ei{AE2).,dG)n\n',\<CN-', 



(5.50) 



e{A{Ei+Ei2),SbG^] 



n\n'g 



< CN' 



We deduce the last inequality by analyzing separately for 51s \ ^'^d 5^2 \ ^^o- 
If 17q C U fii, we can get n = 0. Then 



(5.51) 
(5.52) 



e5y{{AE2)x,G)n2nn[\^Q, 
= 0. 



e{A{E,+E,2),5hG,)n2nn. 

In this case, the modification of (|5.42p is that 

(5.53) e5y{\AS\ + \AE2\,\hG\)r,^^^ 

<Ge5y\\A{S + i?2)||LMr„,..)l|G|lL»(r„.„) 

<Ge^^'^ 5y\\G\\L^(n^\a>^) 

<GN-^ 

where we have used Theorem W/I\ with v ^ I. Collecting (|5.40p . (I5.4ip . (|5.43p . 
(EH, (lOD and we have 

\{eAu,SbG^)\ < C(A^-'' + eiV-^/4)(lniV)|||G'|||. 

If gLflsUfli, using (IOn)) - ([53Dl) leads to 

I {eAu, 6bG,)\ < C{N'P + e'/%){\n N)\\\G\\\. 



□ 
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